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Abstract. We numerically investigate transverse stability and instability of 
so-called cnoidal waves, i.e., periodic traveling wave solutions of the Korteweg- 
de Vries equation, under the time- evolution of the Kadomtsev-Petviashvili 
equation. In particular, we find that in KP-I small amplitude cnoidal waves are 
stable (at least for spatially localized perturbations) and only become unstable 
above a certain threshold. In contrast to that, KP-II is found to be stable for 
all amplitudes, or, equivalently, wave speeds. This is in accordance with recent 
analytical results for solitary waves given in [301 13 1 1 . 



1. Introduction 

The starting point of this investigation is the celebrated Korteweg-de Vries 
equation (KdV) for shallow water waves [26] , i.e. 

(1.1) dtu + ud x u + d xxx u = 0, 

subject to some given initial data u\ t =o = ^o(^) £ ^- Clearly, only accounts 
for uni-directional waves, where in fact one already has taken into account the 
change of coordinate x \-> x + t in order to transform the arising linear transport 
operator 

T±u := d t u + d x u 

into a single partial derivative w.r.t. time, only. The corresponding two-dimensional 
analogue for waves with finite speed of propagation is given by d t u + y/—Au. Here, 
\J — A is the Fourier multiplier corresponding to the symbol ^JQ, + Q , where £ x ,€y £ 

1R denote the frequencies (or Fourier variables) in x and y direction, respectively. A 
generic approach for taking into account weak transverse effects in the propagation 
of quasi uni-directional waves is as follows: By Taylor expanding the dispersion 
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relation 



+ Q « & + ^C 1 ^, when « 1 and 



«1, 



d t u + V— « + -<9 X 1 d yy u 1 



2 ' 

and taking into account the same change of reference frame as before, one formally 
obtains 

1 

2< 

where d~ 1 d yy is defined as a singular Fourier multiplier with symbol — i^ 1 ^. By 
following this approach one formally derives the well known Kadomtsev-Petviashvili 
equation |20j (KP), which after rescaling y ^ y/2y, can be written as: 

(1.2) d t u + ud x u + d xxx u + Xd~ 1 d yy u = 0, t,a;,|/,Gl, A = ±1, 



subject to initial data u\ t =o = U{ n (x,y), such that U[ n (x,0) = uq(x). Equation (1.2) 
can be seen as a (generic) 2 + 1 dimensional generalization of the KdV equation for 
quasi-uni-directional waves with weak transverse effects (see also |15j ). The case 
A = +1 corresponds to the so-called KP-II (or defocusing KP) equation, valid for 
strong surface tension, whereas A = — 1 corresponds to the KP-I (or focusing KP) 
equation, obtained in the case of small surface tension. Two basic conservation 



laws for (1.2) are the mass 



M = u 2 dxdy, 
J JM? 

and the energy 

E = [[ (d x u) 2 -\u* + \(d- 1 d y u) 2 dxdy. 

The KP equation arises, e.g., as an asymptotic model for the propagation of 
capillary-gravity waves [2], or for traveling waves in Bose-Einstein condensates, cf. 
[9]. Analogously to KdV, the KP equation can be written in the form of a Lax 
pair |33j and explicit solutions can be obtained by invoking the inverse scattering 
transform, see, e.g., [H I14j . Generalizations to higher order nonlinearities of the 
form d x u p are also considered (see |32j for a broad review), with p = 4/3 being 
the critical case (in the sense that the potential energy can be controlled by the 
dispersion). Thus, in general, one can expect finite-time blow-up for p > 4/3, cf. 



|27j for more details. Clearly, any ^-independent solution of (1.2) is a solution of 
( |1.1[ ), forming the so-called KdV sector. It is then a natural and physically relevant 
question to ask whether traveling wave solutions of KdV are stable when considered 
as solution of the KP equation. This problem is usually referred to as transverse 
(in-) stability, to be considered in the following for so-called cnoidal waves. 

2. Cnoidal waves 

Periodic traveling wave solutions to the KdV equation can be written in terms 
of elliptic functions. In particular the well-known cnoidal waves are given by 

(2.1) u cn (x, t)=u + I2n 2 k 2 cn 2 (k(x - x - (V + u )t); k), V = An 2 (2k 2 - 1), 

where u^^k^xq g M are arbitrary constants. Here, cn(- ; k) denotes the Jacobi 
elliptic cosine function [7] with elliptic modulus k G [0, 1[. Each solution given by 



(2.1) is periodic in x with period 



P7\ 
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where K(k) is the complete elliptic integral of the first kind. The solution for 
xq = uq = and k = 0.5 can be seen in Fig. [I] 



Figure 1. Cnoidal solution (2.1) to the KdV equation with pa- 



rameters k = 0.5, k = 0.5, and xo = uq = 0. 

Note that in the limit as k —> 1 we have: K{k) — > +oo and cn(x; fc) — )► sech(x) 
in which case the cnoidal solution reduces to the well-known KdV soliton: 

(2.2) Uoo(x, t) = uq + 12ft 2 sech 2 (£ — xo — (4^ 2 + uo)t). 

In contrast to cnoidal waves the KdV soliton is spatially localized, since it expo- 
nentially decays as x — » d=oo. Also note that the amplitude and the wave speed V 
are correlated through the parameter k. In the case of higher order nonlinear it ies, 
similar periodic traveling wave solutions appear, such as dnoidal waves (for mod- 
ified KdV) or snoidal waves. All of them received their name from the fact that 
they can be expressed in terms of Jacobi elliptic functions. 

There has been a large amount of work aimed at understanding the stability of 
traveling wave solutions. The stability of solitary waves under the KdV flow is well 
known since the works of [31 [5]. For periodic solutions, such as cnoidal waves, the 
situation is less complete and there are mainly two types of stability results: Orbital 
(or nonlinear) stability with respect to periodic and sub-harmonic perturbations 
has been shown in pH 110] . and linear and/or spectral stability with respect to 
non-periodic (i.e., localized or bounded) perturbations is proved in [6], 118] . 

This situation becomes even more involved when one asks about the transverse 
stability of traveling waves under the KP flow: Early results in [1 establish linear 
instability of the KdV soliton for KP-I, whereas for KP-II these solitary waves are 
shown to be linearly stable. More recently, rigorous results given in [30, 31 have 
established the existence of a critical threshold in n for KP-I posed on R x x T y , 
where T = R/(27rZ). In particular, for k > k c = 3 -1 / 4 , solitary waves are shown 
to be orbitally unstable for KP-I [30 , whereas for k < n c they are proved to 
be orbitally (nonlinearly) stable |31j . Generalizations of both results to higher 
order nonlinearities are possible. In addition, the full nonlinear stability (including 
asymptotic stability) of KP-II posed on x T y for any kgM has been obtained 
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in |28j . Note that if u(t,x,y) is solution of (1.2), then so is u(—t, —x,y) and hence 
instabilities arising for positive times, also arise for negative times. Numerical 
studies of the stability of the KdV soliton for KP have been performed for instance 
in [24] , 

Concerning the transverse (in-) stability of cnoidal waves under the KP flow, 
several results have been obtained in recent years: Spectral instability for small 
periodic solutions to KP-I posed on x R y has been proved in [XT] , where the 
author allows for two-dimensional perturbations which are either periodic in the 
direction of propagation (with the same period as the one-dimensional traveling 
wave), or non-periodic perturbations. In addition, spectral stability for small solu- 
tions to KP-II on M x x My is also proved in |17] , assuming perturbations which are 
periodic in x and have long wave-length in y. For generalized KP-I type equations 
on R x x Ry, linear instability is proved in [19 . Finally, |16j proves spectral and 
linear instability for KP-I posed on T x x T y , assuming that the perturbation admits 
the same period in x as the cnoidal wave. 

In order to establish a more qualitative picture, we shall in the following perform 
a numerical study of the transverse stability of cnoidal waves under the KP flow. 
The perturbations considered, are either spatially localized in both x and y (see 
Section [4|, or ^/-periodic, and thus non- localized, deformations of cnoidal waves 
(see Section [5|. Before describing our finding, we shall first describe the numerical 
algorithm in the upcoming section. 

3. Numerical Method 



In the present work, the KP equations (1.2) will be numerically studied for 
(x, y,t) GT 2 xR (note that this is the same setting as used in |16j and still close to 
the one of [30, 31 , where y G T.). The periodicity in both spatial coordinates allows 
the use of Fourier spectral methods. The function u(x,y,t) is thereby numerically 



approximated via a truncated Fourier series in both x and y and (1.2) is equivalent 
to 

(3.1) d t u + i^ x u 2 - i^u ± - _^, Q = 0, u\ t=0 = u [n (£ x ,£ y ). 

Here, the singular multiplier — i^ 1 is regularized in standard way as — i(£ x H-zO) -1 , 
see |25j for more details. These methods show excellent convergence properties 
for smooth functions such as the ones studied in here. In addition, they have the 
advantage that only minimal numerical dissipation is introduced, which is impor- 
tant since numerical dissipation could suppress some of the dispersive effects we 



are interested in. In view of (3.1), we see that the use of Fourier methods forces 
us to solve a large system of ordinary differential equations in time. Due to the 
third derivative in x and the non-local term this system is in general very stiff. 
However, since the stiffness is in the linear part only, there are efficient methods to 
allow for time integration of high precision which avoids a pollution of the Fourier 
coefficients with numerical errors. In [22] and |23j it was shown that fourth order 
exponential integrators are very efficient for KdV-like equations and we shall use a 
method of Cox and Matthews [8] in the following. 

In the following, we shall always use N x = 2 10 modes in x-direction and N y = 
2 8 modes in ^/-direction for x G [— ttL x , ttL x ] and y G [— 7rL y , 7tL y ]. We choose 
L x = 87rK(k) I 'ft, where K(k) is again the complete elliptic integral of the first 
kind. This choice ensures periodicity of the solution. For the ^-direction we choose 
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L y = 2. In all our numerical simulation, we set in (2.1) uo — xo = and k — 0.5 
so that the only remaining parameter is k £ R. The elliptic functions and integrals 
are computed with Matlab commands to machine precision. 

To test the quality of the numerical code, we first propagate the cnoidal solu- 



tion (2.1) with the full KP code. As mentioned before, the KP equation admits 
mass conservation, i.e., conservation of the L 2 norm, but since this feature is not 
implemented in the code, numerical errors will always lead to a time-dependent 
L 2 -norm. The latter can be used as an indicator for the quality of the numerical 
results, see also the discussion in [23 . As a first test case, we first consider the 
KP-I equation with initial data 

u in (x,y) = u cn \ t =o(x), and k = 2. 

We find that the L 00 -norm of the difference between numerical and exact solution 
is of the order of 2.6 x 10 -8 for N t = 10 4 time steps. In addition the numerical 
error in the conservation of mass 

A-=l-^ 

L 2 (0)' 

where L 2 (t) is the numerically computed (discrete) L 2 norm at time t £ M, is 
roughly 10 -10 . For KP-II, with the same initial data, the values are found to be 
almost identical. This shows that the quantity A overestimates numerical precision 
as discussed in |23j by roughly two orders of magnitude and can be used as a 
reliable indicator of numerical accuracy. We shall always aim for values of A of the 
order of 10 -5 or smaller in order to ensure an accuracy much better than plotting 
accuracy. For smaller amplitudes, i.e., k < 2, the accuracy is found to be even 
higher: For k = 0.5, we find A = 10 -12,6 and for the L°°-norm of the difference 
between numerical and exact solution 1.87 x 10 -13 for KP-I, and A = 10 -12,6 and 
1.4 x 10 -13 for KP-II. Working with double precision in Matlab, which allows for 
a precision of 10 -14 , we conclude that the solution can indeed be propagated with 
machine precision. 

4. Gaussian perturbations 



In this section we shall consider the KP equation (1.2) subject to initial data 
of the form 

(4.1) u in (x, y) = u cn \ t =o(x) + Up(x, y), 

where u cn is the cnoidal wave for t = (and uo = xo = 0) and u p is a (small) 
perturbation which is assumed to be localized with respect to both x and y. 

Next, we recall the well known fact (see for instance [29j for a rigorous proof) 
that KP equations satisfy the following constraint, in the sense of Riemann Inte- 
grals: 



(4.2) / u(t,x,y)dx = Q, y eX,t^0, 

Jx 



where X = R, or X = T, respectively. Equation (4.2) holds, even if the initial 
condition U[ n does not satisfy the constraint, see |25j for a numerical example. In 
such a situation the solution ceases to be continuous at t = 0, however. Since this 
leads to numerical problems, we will only consider initial conditions which satisfy 
the above constraint. This is obviously the case for the cnoidal KdV solution, but 
we shall also impose it on the perturbations u p . A possible way to do so is to 
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consider U[ n = d x f(x,y), where / is a periodic or Schwarzian function, which can 
be considered as essentially periodic if the period is chosen large enough such that 
the solution at the limits of the computational domain decreases below machine 
precision. The latter is important to avoid the appearance of Gibbs phenomena due 
to discontinuities. 

To be more precise, we shall study localized perturbations of the following form 



(4.3) 



Up(x,y) =xe ^ 2 +^ 2 ), 



i.e., an x-derivative of a Gaussian. With the same parameters as above and k = 2, 
we obtain for the time evolution of the initial data ( |4.1| ) that at t — 2 the relative 
L 2 norm is A ~ 10 -10 . It can be seen in Fig. [2] that the solution stays close to 
the cnoidal form for some time, but will finally disintegrate into a periodic array of 
humps. It is conjectured that for (x,y) G M 2 , solutions to KP-I asymptotically de- 



t=o.o 



t=0.7 





t=1.4 



t=2 





Figure 2. Solution to the KP-I equation (1.2) for several values 



of t obtained for initial data given by a cnoidal wave (2.1) with 



parameters k = 2, k = 0.5, and xq — uq =0, and a Gaussian 



perturbation u p (x,y) of the form (4.3) 



compose into several lump solitons, i.e., solitons which are localized in both spatial 
directions with an algebraic decay towards infinity. Numerically, this was already 
observed in |24j . However, there is no periodic analogue of the lump soliton known, 
and thus it is not clear whether there are doubly periodic solutions to KP-I which 
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appear as asymptotic solutions of a perturbed cnoidal solution. There are, how- 
ever, doubly periodic solutions to the KP equations, which can be given in terms 
of multi-dimensional theta functions on hyperelliptic Riemann surfaces, see, e.g., 
113] . The simplest such example is a travelling wave on a genus 2 surface that, 
which for KP-II is depicted in Fig. [3] 




Figure 3. Doubly periodic hyperelliptic KP-II solution of genus 
2 for a fixed time. The profile travels with constant speed in re- 
direction. 



The situation is considerably different if the same initial data ( |4.3[ ) are used 
for the KP-II equation (for which we essentially obtain the same value of A). Note 
that even though the equation is defocusing, the perturbation cannot be radiated 
away to infinity as would be the case for the perturbations of the KdV soliton on all 
of R 2 . In our periodic setting, they can just smooth out over the whole period or 
distribute periodically. In Fig. [4] w e show the difference between the KP-II solution 
corresponding to initial data (4.1) with u p = (i.e. a purely cnoidal solution) and 
u p given by (4.3), respectively. It can be seen that, as time increases, the initial 
hump in fact smoothes itself out over the whole domain. 

So far, all simulations were done with k = 2 which is above the critical value 
k c established in [31 for the transverse stability of solitary wave (obtained for 
k — » 1 in the cnoidal solution). Our simulations confirm that for n > n c the cnoidal 
solution is unstable for KP-I, but stable for KP-II. To check whether there is a 
similar threshold below which cnoidal waves become stable for KP-I, we shall also 
consider the value k = 0.5 which is well below k c . To compare similar situations, 
we multiply the final time by a factor of 16 (since the the velo city V in (2.1 ) scales 
as n 2 ) and in addition consider a perturbation of the form (4.3) but rescaled by 
a factor u p /16. The latter ensures that the perturbation is of the same relative 
amplitude as in the case k = 2 considered above. The code produces for both KP-I 
and KP-II an error A < 10 -12 . In Fig. [H] it can be seen that the KP-II situation is 
essentially unchanged in this case, i.e., the cnoidal solution is stable for KP-II also 

for k < Kn. 
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Figure 4. Difference of the solution to the KP-II equation (1.2) 
for the initial data (4.1) with u p given by(4.3), respectively. The 
cnoidal initial data are chosen with n = 2 > k c , k = 0.5 and 
x =u = 0. 



The corresponding situation for KP-I is found to be considerably different when 
compared to the situation shown in Fig. [2] Indeed, it can be seen in Fig. [6] that for 
k = 0.5 < n c the solution is almost identical to the KP-II case, and thus also stable. 
It consequently appears that, just as for solitary waves, there is also a critical value 
for n separating the stable from the unstable regime in KP-I. 



5. Periodic deformations 



In this section we shall consider initial perturbations of (2.1), which are no 
longer spatially localized. Rather, we shall consider a ^/-periodic deformation of the 
of the form 

(5.1) u- m (x,y) = 12k; 2 /c 2 cn 2 (x + Scos^J ]kj , 5 > 0., 

For (5 = this yields the usual cnoidal wave (with xq = uq = 0). In the following, 
though, we shall choose S = 0.4 yielding a strongly deformed cnoidal type solution 
as initial data. Note that these kind of ^/-oscillatory perturbations are different 
from the x-periodic perturbations considered in, e.g., \16\ 117] . 
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Figure 5. Difference of the solution to the KP -II e quation (1.2) 
for initial data u- in = u cn + jqU p with u p given by (4.3), respectively. 
The cnoidal initial data are chosen with n = 0.5 < ft c , k — 0.5 and 
= = 0. 



For small n = 0.5 these initial data lead a solution of KP-II equation as shown 
in Fig. [7| In fact the solution oscillates around the cnoidal solution which thus 
appears stable in this sense. Indeed, these oscillations can be seen best by tracing 
the L°° norm of the solution as a function of time, see Fig. [5] One can expect these 
oscillations to go on indefinitely. For the KP-I, the same initial data lead to an 
essentially identical behavior as can be seen in Fig. [8] and Fig. [9] Thus for small k, 
KP-I and KP-II show very similar behavior for perturbed cnoidal initial data. 

Analogous to the case of localized perturbations, the situation becomes different 
for larger For k = 2, the KP-I solution corresponding to initial data (5.1) is 



given in Fig. 10 The solution appears to be periodic in time in the sense that 
there are times where it is close to the initial configuration. But there is much 
more structure in this type of doubly periodic solution than in the genus 2 example 
shown in Fig. [3] In particular there seems to be a 'breathing' structure similar to 
the so-called breathers in the nonlinear Schrodinger equation (NLS), see, e.g., |12j 
(for the relation between NLS and KP solutions see O I21j). 

For KP-II the situation for n = 2 is similar to the one of KP-I as can be 



seen in Fig. 11 Again, the initial configuration reappears from time to time, but 
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Figure 6. Difference of the solution to the KP-I equ ation ( 1.2 ) for 
initial data U[ n = u cn + -^u n with u n given by (4.3), respectively. 

0.5 < k c , k = 0.5 and 



^cn n" i6 vviLja u, p 5iv^ 

The cnoidal initial data are chosen with n 



x =u = 0. 



there are very rich oscillations patterns. It thus seems to be that for such kind 
of de-localized (^/-periodic) perturbations, both KP-I and KP-II exhibit a similar 
behavior and both are not stable in the usual sense. 



6. Conclusion 

We presented numerical simulations concerning the transverse stability and 
instability of cnoidal waves for KP equations. By means of a spectral method we 
found that: 

• For spatially localized perturbations KP-II is stable for general k G R, 
whereas KP-I is only stable for k below a certain threshold (and unstable 
above). This is in agreement with recent rigorous results on the transverse 
stability of solitary waves (obtained from the cnoidal solution in the limit 

• For ^/-periodic perturbations and small ft, both, KP-I and KP-II are stable. 
For large k both KP-I and KP-II appear to be unstable in the usual 
sense, as both exhibit a doubly periodic type solution with rich oscillatory 
behavior. It remains a challenging open problem to give a more precise 
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t=0 



t=2.8 





t=5.6 




V " 5 -20 




V " 5 ^ -20 



Figure 7. Solution to the KP-II equation for several values of £, 



obtained for initial data given by (5.1), i.e. periodically deformed 



cnoidal initial data with n = 0.5, k = 0.5 and xq = uq = 0. 



f I i 



f\(\t\ /i (In /i d/1 



Figure 8. L°° norm of the KP-II solution shown in Fig. [7] as a 
function on time on the left and for KP-I of Fig. [9] on the right. 

description of the long time behavior of cnoidal waves under such type of 
deformations. 
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